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We show that a one dimensional device supporting a pair of Majorana bound states at its ends 
can produce remarkable Hanbury Brown-Twiss like interference effects between well separated Dirac 
fermions of pertinent energies. We find that the simultaneous scattering of two incoming electrons or 
two incoming holes from the Majorana bound states leads exclusively to an electron-hole final state. 
This "anti-bunching" in electron-hole internal pseudospin space can be detected through current- 
current correlations. Further, we show that, by scattering appropriate spin polarized electrons from 
the Majorana bound states, one can engineer a non-local entangler of electronic spins for quantum 
information applications. Both the above phenomena should be observable in diverse physical 
systems enabling to detect the presence of low energy Majorana modes. 



I. INTRODUCTION 

Quantum Indistinguishability has striking manifesta- 
tions when two identical particles are brought together 
at a beam splitter. For example, two bosons in identi- 
cal states would "bunch" together when exiting a beam- 
splitter purely due to interference effects [1]. Two 
fermions, on the other hand, would exit separately or 
"anti-bunch" These effects are indeed an instance of 
the celebrated Hanbury Brown-Twiss effect, which has 
recently also been tested with Helium atoms [||. The 
same quantum indistinguishability is exploited for the 
production of entangled photons Q, and can also be 
used to entangle generic massive particles Q. Of course, 
all these effects can occur only when the particles are 
brought together spatially, for instance, at a beam split- 
ter. It is thereby interesting to look for settings where 
rather well separated identical particles could manifest 
such phenomena. 

Here we report on the possibility of engineering a non- 
local beam splitter enabling the above class of phenomena 
for distant charged fermions. Here, by "non-local" we 
mean spatially extended. Going beyond the usual two 
particle interference in orbital/momentum space, here 
one finds a Hanbury Brown-Twiss effect in the electron- 
hole internal pseudospin space. This is enabled by Majo- 
rana mid-gap low energy modes which transform between 
electrons and holes [6|, effectively making them indis- 
tinguishable in a scattering experiment. This Hanbury 
Brown-Twiss effect is thereby a detector of the Majorana 
modes. 

Recently, low energy Majorana (neutral charge self- 
conjugated fermion) modes located at the edges of linear 
devices have been predicted to induce non-local phenom- 
ena 7 9]. Indeed there are a variety of platforms to re- 
alize such devices such as a quantum wire immersed in a 
p-wave superconductor 0, , cold-atomic systems mim- 



icking p-wave superconductors |lT| , to pol ogical insulator- 
superconductor-magnet structures |§, ll2l. [l3| and poten- 
tially also semiconductor systems [lj, [l5[ . The evidences 
of their non-local nature are distance independent tunnel- 
ing Q, crossed Andreev reflection H and teleportation- 
like coherent transfer of a fermion [9j . Finally, they may 
be easily manipulated @ and are relevant excitations also 
in conventional superconductors [l6j . So far, the primary 
application envisaged for these fermions has been topo- 
logical quantum computation pj}. As the second key 
result of this paper we will show another use of these 
modes, namely that Majorana bound-states (MBS) could 
be used to engineer entanglement between the spins of 
well separated particles, a pivotal resource in quantum 
information. 

The paper is organized as follows. In Section II we con- 
sider the scattering of Dirac fermions off the edge MBS of 
energy Em in the spinless model investigated in Q. Here, 
we show that, when the energy of the incoming fermions 
is nearly resonant with Em, the edge MBS induce a beam 
splitting process which acts like an equally weighted 4- 
port beam splitter, with ports corresponding to both spa- 
tial and electron-hole isospin states. In Section III we 
show that with two incident Dirac fermions, this allow 
for fermion antibunching in the pseudospin space which 
has holes and electrons as its two states. This is one of the 
central results of our paper. In Section IV we determine 
the signature of this fermion antibunching in the zero- 
frequency spectral density of the current fluctuations in 
the leads. Section V generalizes the results of section II 
by accounting for the spins of the fermions. In section 
VI, we show that the edge MBS allow for generating en- 
tanglement between the spins of distant electrons only 
by pertinently choosing the polarizations of the incom- 
ing fermions. In Section VII we analyze a few condensed 
matter settings where our findings may be helpful in de- 
tecting the presence of MBS; section VIII is devoted to a 
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since- for the evaluation of the S-matrix- the leads my be 
effectively removed by introducing complex embed ding 
potentials 18]. Using the approach developed in fl8Lll9f 
one may then write the unitary scattering matrix S in a 
form which is formally independent on the model used to 
describe the leads; for the scattering of fcrmions from the 
MBS located at the the ends of the quantum wire shown 
in Fig[T] the S matrix has been computed in Q as 

S(E) = 1 + 2mW\H M -E- mWW^W (2) 



where W is a rectangular matrix 



FIG. 1: Non-local beam splitter and electron spin entangler. 
The MBS are shown as empty ellipses 1 and 2. One specific 
realization where MBS occur at boundaries between magnets 
(M) and superconductors (SC) deposited on quantum spin- 
Hall insulators is depicted, though our results hold more gen- 
erally. Incoming and outgoing particles are shown by arrows, 
and may, in practice, be tunneled in/out by STM tips or elec- 
tron pumps acting as leads. 



few remarks on our results. For a reader interested more 
in the logical steps leading to our central results, rather 
than the full technical details, we recommend focussing 
primarily on sections III and VI, taking the relevant scat- 
tering matrices from the sections immediately preceding 
them. 



II. NEARLY RESONANT ELECTRON (HOLE) 
SCATTERING FROM MAJORANA EDGE 
STATES 

We consider a one dimensional device supporting two 
weakly coupled MBS at its ends as shown in FigQ] The 
MBS are labeled as 1 and 2 and schematically shown as 
empty ellipses in the figure. As the separation between 
the MBS increases their energy Em decreases exponen- 
tially [7j. For the sake of clarity, we will first show how 
this device produces Hanbury Brown- Twiss like interfer- 
ence effects between spatially separated Dirac fermions in 
the spinless models investigated in 0, [Hj; later, we show 
how all the results are valid for more realistic spinfull 
physical settings [HI, [3 EH . The Hamiltonian describ- 
ing the weak coupling between the MBS 71 and 72 is 
given by 

Hm = iE M lil2, (1) 

where jj are Majorana operators defined by jj = 7! 
and satisfying jjjk + 7fc7j = "28kj (in our definition, the 
jj = Cj + Cj in terms of Dirac fermionic operators Cj). 

Leads, also labeled as 1 and 2, are connected to the de- 
vice as shown in Fig|TJ allowing for the scattering of Dirac 
fermions (electrons or holes) from each of the MBS; we 
further assume that lead 1 is coupled only to bound state 
1 while lead 2 only to bound state 2. The Hamiltonian 
describing the leads needs not to be specified at this stage 
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in the basis {|ei), \e2), \h%), |/i2)}> with \ej) and \hj) rep- 
resenting an electron and a hole in the lead j. W de- 
scribes the coupling of the scatterer (Hm) to the leads 
and E is the energy of the incident electrons/holes. The 
entries wj of the W matrix are related to the couplings 
to the leads Tj = 2-kw] 0. 

For our purposes, it is convenient to assume that E » 
Tj, as well as E w Em (i.e., the energies of the incoming 
Dirac fermions are tuned to be nearly resonant with the 
Majorana coupling energy). Under these circumstances, 
the S matrix simplifies to 



S 




(3) 



where the basis is, again, {|ei), |e 2 ), \h 2 )}- Note 
that this regime is different from the one considered 
by Akhmerov et. al. [8j, where only the terms corre- 
sponding to crossed Andreev reflection (i.e (ti2\S\ei) and 
(hi\S\e2)) are maximized. Here, we work in a regime 
where all the entries of S have the same magnitude. It 
is the implications of this scattering matrix S of Eq.Q 
that we work out in this paper. It is this S that enables 
both the non-local Hanbury-Brown Twiss interferometry 
in isospin space and the non-local entanglement genera- 
tion. 

Let us first illustrate the action of the above S ma- 
trix by describing what happens to electrons tunneling 
in from one end of the one dimensional device. If, at 
time t, a single electron tunnels into the Majorana mode 
located at site 1, i.e., the incoming state is c||0), it trans- 
forms, under S, to 
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where c 1 



j (dj) creates an electron (hole) at site j. In 
Eq.dU, we have used MBS above the arrow to indicate 
that Majorana bound states are responsible for the pro- 
cess. Since the transformation ^ is equivalent to a four 
port beam splitter, with MBS inducing the beam split- 
ting process, one can equally well take MBS to stand 
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for "Majorana Beam-Splitter". Eq.© implies that an 
incoming electron has ^th probability of coming out of 
each site as an electron or a hole. If another electron 
scatters at a different time t' on the Majorana mode lo- 
cated at position 2, it will also scatter with exactly the 
same probabilities for the four possible outcomes. The 
joint probability for two incoming electrons to exit as 
two electrons or two holes (whichever the output port) 
would thus be i. 



III. HANBURY-BROWN EFFECT IN 
PSEUDO-SPIN SPACE 

We will now show that when t = t' , i.e., simultane- 
ous scattering, two particle interference can take place 
so that the probability of two electrons or two holes ex- 
iting is completely suppressed. By t — t' we mean that 
the wavepackets of the two incoming electrons (holes) are 
large enough so that their time of arrival cannot be dis- 
tinguished when one observes them after the scattering. 

When two electrons scatter simultaneously one at site 
1 and the other at site 2, one has 
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= ^{ic\d\ - c[dl + c\d\ + icldl)\0). (5) 

In the last step of Eq.©, we have used d^(E) — Cj(—E) 
(which effectively embodies the indistinguishability be- 
tween an electron and a hole), where E is energy. From 
Eq.© one sees that the probability for two outgoing 
electrons (holes) after the scattering is zero. Exactly 
the same holds when two holes scatter simultaneously at 
leads 1 and 2. This is an interference effect in the same 
sense as the anti-bunching of fermions at a normal two 
port beam splitter, where fermions cannot exit through 
the same port. Instead of being in the spatial channels, 
here the anti-bunching is in the internal pseudospin space 
which has particle and hole as its two states. The uni- 
tary conversion of an electron to a hole, is, per se, not 
surprising in view of Refs.Q. 

Of course, in a practical realization, the condition 
E ~ Em required for obtaining the scattering matrix 
S of Eq. (|3j) may not be exactly met. To see the effect 
of an energy mismatch, we denote by 5E the amount by 
which E deviates (either positively or negatively) from 
Em] this deviation is, however, assumed to be much lower 
than Em itself (i.e., SE « Em)- Without assuming 
SE << Em, one may end up in qualitatively different 
regimes: e.g., for SE comparable to —Em, one reaches 
the regime of Ref.[8| of only crossed transmission. For 
SE << Em, the scattering matrix as a function of SE is 
given by 

ci *T SE 

= jETiT 6~E~+iT (6) 



where r = Tx ~ T2 and / is the 4x4 identity matrix. 
In deriving Eq.©, one ignores the second and higher 
powers of both SE/E M and T/E M as E M » AE,T. It 
is easy to check that, despite the above approximation, 
Sse is unitary; furthermore, Eq.© holds for any value 
of the ratio SE/T as long as Em » AE,T. Using Sse, 
one readily obtains that the probability of observing an 
electron-electron output state becomes finite and equal 



to 



(SEV 



which, of course, vanishes when E ~ E 



M ■ 



(SE) 2 +T 2 

Before ending this section, as a brief aside, we point 
out that, when one electron and one hole scatter at sites 
1 and 2 respectively, for SE « T, the incoming state 
44|0) evolves to |(— c\c\ — ic\d\ +ic\d\ — d\d\)\0), im- 
plying the interferometric vanishing of the probability of 
one outgoing electron and one outgoing hole in separate 
leads. We will not delve further into this case, but next 
proceed to discuss the signatures of the Hanbury-Brown 
interferometry in the case of two incident electrons. 



IV. SPECTRAL DENSITY OF CURRENT 
FLUCTUATIONS: A SIGNATURE OF FERMION 
ANTIBUNCHING IN PSEUDOSPIN SPACE 

In the previous section we have described the scatter- 
ing as a process where one sends particles one by one 
through the leads at specific times. However, in practice, 
rather than controlling times, one could control the ener- 
gies ei and £2 of the particles in their respective leads, so 
as to make them behave indistinguishably when e\ ~ £2. 
Then, the standard way to observe the predicted fermion 
anti-bunching is through a measure of the correlations 
between the currents in leads 1 and 2. The current in 
lead j may be written as [2(], Hl[ 
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where aj and bj denote the incoming and outgoing par- 
ticles and £ and £ are the energies of the particles and 
v is the density of states of the incoming electrons. The 
spectral density of the current fluctuations 81 j = Ij — (Ij ) 
between the leads at zero frequency is [2l[ 

hp f T 

Pij = lim T ->oo -=■ / dt Re (Sh (t)SI 2 (0)) . 
1 Jo 



(8) 



Using Sse of Eq.© and considering an incoming two 
electron state c|(ei)<4(e2)|0), where ct(ej) denotes an 
electron of energy ej in lead j, one finds 



hv{(SE) 2 +r 2 } ; 



: {(SE) 2 -T 2 }S eu 



(9) 



where S ei>e2 is the Kronecker delta function. Note that, 
when the incident electrons are distinguishable i.e., £1 ^ 
£ 2 , then, as expected, — since for an electron 
exiting one lead there could equally well be an elec- 
tron or a hole exiting the other lead. When, instead, 
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ei = E2 (i.e., the particles are indistinguishable), then for 
\SE\ < |r|, the domination of the electron- hole final state 
(as in Eq.©) makes P^j < 0, which allows to the detect 
the predicted "anti-bunching" in pseudospin space. For 
\SE\ > |r| a process of amplitude TSE in which only one 
of the electrons scatter, while the other remains in its 
lead, dominates; Fermi statistics now makes the electrons 
anti-bunch spatially (the more conventional antibunch- 
ing [1, HU), contributing to a positive Py. As in Ref.jsJ], 
our results are not inconsistent with those of Bolech and 
Demler [22| , since their results apply when the energy of 
the incoming electrons is much higher than Em- 



V. SCATTERING MATRIX IN THE SPINFULL 
CASE 

So far our analysis has been confined to the spin- 
less model investigated in 0, while for the promis- 
ing implementations [l2T[T5| , the Majorana modes should 
involve superpositions of operators of different spins. 
For example, for a realization in a ferromagnet-s-wave 
supcrconductor-ferromagnet structure on a quantum 
spin-Hall edge [l3j . one has 
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(10) 



where Cj_ a creates an electron with spin 
a in lead j. Defining the spin states 

| ± y) = -j=(\ f) i i\ I-))) an d using the basis 

{|ei,+ v ), \e-2-y), \hi, +y ), \h 2 ,- y ), \ei.-y), \e2,+ v ), \hi,- y ), 
\Ji2,+y)}, the scattering matrix is found to be 
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where in (TTTI) . / and are the 4x4 Identity and null 
matrices, while S is the scattering matrix given by Eq. (|3]) . 

When one uses S^nfuii to study the scattering of the 
incident state c\ _ y 4 +J) |0), one only needs the lower- 
right 4x4 block of Sgpinfuii. Thus, precisely the same 
electron- hole output state as in Eq.© is obtained, apart 
from the fact that, now, the spin indices —y and +y 
are pinned to the sites 1 and 2 respectively. Thus, by 
choosing the spin polarizations of the incoming electrons 
pertinently, one can observe all the effects described 
till now. This should be possible in a variety of sys- 
tems as Majorana modes of the form given by Eq. (fT0)) 
are quite generic, e.g., also realizable in semiconductor- 
superconductor- magnet structures [l5j . 



VI. ENTANGLEMENT OF DISTANT 
ELECTRON SPINS 

We now propose a protocol for the generation of entan- 
glement between spins of well separated particles incom- 



ing at site 1 and at site 2. For this purpose, we choose 
the realization of Majorana fermions given by Eq. fTU]) 
and make two electrons with parallel spins in the f direc- 
tion come in simultaneously i.e., choose the initial state 



'l,f 2,t 



d, JO). Then, using S sp i n f u u, one gets 
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as 
which 



lm- 
are ob- 
state is 



where ... denotes terms such 
c] c] ,,ctcl ,,c\„d] , and d\ d\ ,, 

are not relevant to our discussion. Eq, 
plies that, when two outgoing electrons 
tained in leads 1 and 2 separately, their 
I0u = 7g(l t)i| t)a - U)i| 1)2 + 2| t)i| ;)a) where, as 
it is usually done d, [2lJ , one uses the lead label to label 
the electron. \£)i2 is an entangled state of the spins 
of electrons 1 and 2, with the amount of entanglement 
(as quantified by the von Neumann entropy of one 
of the particles [Hj]) being 0.19 ebits. Though the 
entanglement is not very high, |£}i2 is a pure state, 
and hence of value in quantum information, as its 
entanglement can be concentrated without loss by local 
means [241] . Moreover, the probability of obtaining two 
outgoing electrons in separate leads (i.e., |£)i2) is rather 
high, namely 3/8. At the expense of decreasing this 
probability, one may improve the degree of entanglement 
of the generated state by tuning the polarizations of the 
incoming electrons. For instance, if the incoming state 

is (vTo c l + y + 7Io c U)(^o4, +! , + 7b4,- H )|0), one 
obtains an output state of entanglement 0.75 ebits, while 
the probability of the generation this state becomes 
0.055. The spin entanglement of the outgoing electrons 
could be measured by passing them through separate 
spin filters as in Ref. [23l ]. 

Unlike the entanglement generation scheme of Ref. Q , 
here particles polarized parallel to each other suffice to 
generate entanglement. Importantly, in our protocol, 
particles at a distance from each other can be made en- 
tangled; this may avoid the decoherence arising necessar- 
ily from the transport needed to separate the particles 
after a local entangling mechanism. In addition, the dis- 
tance between the entangled particles can be enhanced 
by putting n copies of our setup in series with leads con- 
necting the end of one copy to the beginning of another. 
The probability of obtaining the state |£)i2 in the left- 
most and rightmost leads will then be (3/8) n . 

One might think that in analogy with Ref. 21], per- 
haps it is possible also to detect entanglement between 
distant electronic spins by injecting to the opposite ends 
of our one dimensional device. As a maximally entan- 
gled state of two spins can be written in any basis, let 
us consider the | ± y) basis for spins. In this basis, 
two of the maximally entangled states can be written 

as |V' ± )i2 = (4,+i,4 ,-y ± c l -y c 2,+ y )\ )- The detection 
of entanglement at a normal 50-50 beam splitter relies 
crucially on both the incoming states c\ +y c\ _ |0) and 
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cj + |0) evolving at the beam-splitter and interfer- 
ence (i.e., cancellation/addition) between the terms re- 
sulting from the evolution each of the above two states. 
Only as a result of these cancellations/additions does the 
bunching/ anti-bunching effects evidencing entanglement 
arise. However, here the term c\ + c| _J0) does not even 
evolve under the action of S^infull so that interferences 
are impossible. Thus though our device can generate spin 
entanglement it cannot detect spin entanglement. 



VII. A FEW CONDENSED MATTER SETTINGS 

One simple setting where the non-local two particle 
interferometry and the entanglement generation between 
distant electrons from MBS may be observed can be en- 
gineered with magnet-superconductor-magnet junctions 
deposited on the edge of a 2D quantum spin Hall insula- 
tor [l,[l3j]. Just as in Ref.Q, one can observe these effects 
when the Majorana modes are separated by a distance 
d of several micrometers at temperatures of the order of 
10 mK. For this setting, the explicit form of the Ma- 
jorana operators is exactly the same as in Eq. ([TU|) [l3j . 
Interestingly, strong spin-orbit coupled quantum wires in 
proximity with ferromagnets and superconductors also 



_given in Eq. flOj) 

22], also in these 



support the realization of MBS 1 
|15| . As in previous proposals [131 [1 
settings, two STM tips could act as the leads 1 and 2 
to observe the non-local two particle Hanbury Brown- 
Twiss interferometry. For the entanglement generation, 
instead, it will be more useful to have synchronized elec- 



tron pumps [25[ feeding in the incoming electrons. In 
addition, the filtering of the desired state |£)i 2 can be 
achieved by pumps capturing exactly one outgoing elec- 
tron from each Majorana bound state. 



VIII. CONCLUSIONS 

In this paper, we showed that a one dimensional de- 
vice with two Majorana bound states at its ends yields 
a Hanbury-Brown Twiss effect in the internal electron- 
hole pseudospin space which may be detected in real- 
istic condensed matter settings through current-current 
correlations. This is a departure from all the known 
multi-particle interference effects which have manifested 
themselves in spatial bunching and antibunching or spin- 
spin correlations. Fundamentally, it can be regarded as 
a manifestation of the quantum indistinguishability be- 
tween electronic annihilation and hole creation evidenc- 
ing the presence of Majorana bound states. The same 
settings may also be used to engineer a non-local entan- 
gler of distant electronic spins, which may enable cir- 
cumventing the decoherence arising from the transport 
needed to separate entangled particles. 

Acknowledgements: We thank C. Chamon, R. Egger, 
R. Jackiw, M. Pepper, S. Y. Pi, G. W. Semenoff and S. 
Tewari for fruitful discussions. SB (PS) thanks the Uni- 
versity of Perugia (University College London) for hospi- 
tality and partial support. SB thanks the EPSRC, UK, 
the Royal Society and the Wolfson Foundation. 



[1] C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. Rev. Lett. 

59, 2044 (1987). 
[2] R. C. Liu et. al, Nature 391, 263-265 (1998). 
[3] T. Jeltes et. al., Nature 445, 402 (2007). 
[4] Y. H. Shih and C. O. Alley, Phys. Rev. Lett. 61, 29212924 

(1988). 

[5] S. Bose and D. Home, Phys. Rev. Lett. 88, 050401 
(2002). 

[6] L. Fu and C. L. Kane, Phys. Rev. Lett. 102, 216403 
(2009); A. R. Akhmerov, J. Nilsson and C. W. J. 
Beenakker, Phys. Rev. Lett. 102, 216404 (2009) 
[7] G. W. Semenoff and P. Sodano, J. Phys. B: At. Mol. Opt. 

Phys. 40, 1479 (2007). 
[8] J. Nilsson, AR. Akhmerov, and C.W.J. Beenakker, Phys. 

Rev. Lett. 101, 120403 (2008). 
[9] L. Fu, Phys. Rev. Lett. 104, 056402 (2010). 
[10] A. Yu. Kitaev, Phys.-Usp. 44, 131 (2001). 
[11] S. Tewari et. al, Phys. Rev. Lett. 98, 010506 (2007). 
[12] L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 0964407 
(2008). 

[13] V. Shivamoggi, G. Refael, and J. E. Moore, Phys. Rev. 

B 82, 041405 (2 010). 

[14] J. D. Sau et. al, arXiv:1006.2"829V 2 (2010); J. D. Sau et. 



al, Phys. Rev. Lett. 104, 040502 (2010). 

[15] Y. Oreg, G. Refael, F. von Oppen, larXiv:1003.PT45V 2 
(2010). 

[16] C. Chamon et. al, Phys. Rev. B 81, 224515 (2010). 
[17] A. Yu. Kitaev Ann. Phys. (NY.) 303, 1 (2003); C. Nayak 

et. al, Rev. Mod. Phys. 80, 1083 (2008). 
[18] S. Datta Electronic Transport in Mesoscopic Systems, 

(Cambridge University Press, 1995), pp. 145-157; D.S. 

Fisher and P. A. Lee, Phys. Rev. B 23, 6851 (1981). 
[19] T. Ghur, A. Miiller-Groeling, H. Weidenmuller, Phys. 

Rep. 299, 189-425 (1998). 
[20] M. Biittiker, Phys. Rev. Lett. 65, 2901 (1990). 
[21] G. Burkard, D. Loss, E. V. Sukhorukov, Phys. Rev. B 

61, R16303 (2000). 
[22] C. J. Bolech and E. Demler, Phys. Rev. Lett. 98, 237002 

(2007). 

[23] O. Sauret, T. Martin, D. Feinberg, Phys. Rev. B 72, 

024544 (2005). 
[24] C. H. Bennett et. al, Phys. Rev. A 53, 2046 (1996). 
[25] M.D. Blumenthal al, Nat. Phys. 3 343 (2007); S. J. 

Wright et. al, Phys. Rev. B 80, 113303 (2009). 



